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INTRODUCTION 

Scarce resources and competition lead to deeper improve- 
ments in the productivity, effectiveness, and efficiency of 
all the major activities in production/manufacturing. Some 
improvements are obtained through better product design, 
materials, equipment, and instrumentation (sensors and 
controls), in the realm of engineering; and some by making 
best use of the assets that already exist, partially achieved 
by experience and trial-and-error, but truly completed 
in the realm of technical management, process model- 
ing, and operations research (as of the field in industrial 
engineering). 

There are many activities or domains in production/ 
manufacturing that are practical and worthwhile to improve 
and optimize, for example, product design, process design, 
purchasing, regulatory control, supervisory control, mainte- 
nance, plant engineering, scheduling, production planning, 
supply chain, shipping, storage, and several others. Each of 
these areas requires a particular approach to improvements 
and optimization and its associated technology. What is prac- 
ticed, of course, is practical optimization, balancing extra 
efforts with marginal benefits obtained. 

This chapter focuses on operation optimization during 
regular production with the existing facilities, in particu- 
lar, on the supervisory control decisions of process-setting 


adjustments, which are the directions given to the process 
control system to perform its regulatory and more basic 
supervisory control. The challenge in production manage- 
ment is how to best operate the existing assets, that is, not to 
waste process capabilities, as aspect of excellence and lean, 
to maximize operating performance. 

Operating performance includes first and foremost reli- 
ability meeting minimum requirements such as safety, 
protecting process life, regulations, and commitments to 
customers and the firm (quality specs, volume, and budget). 
For many production operations, achieving this reliability 
is deemed a sufficient success. However, there are many 
production managers and “cultures” who also strive to con- 
tinually expand the firm’s competitiveness and readiness 
for future demands, such as, by producing more consistent 
quality around nominal blueprint/formula values (to satisfy 
customers better and thus increase future demand), increas- 
ing production volume (when sales can support it, one of the 
greatest contributions to extra profitability that manufactur- 
ing can deliver), reducing costs and losses below budget, and 
reducing emissions below government regulations. These 
aspects of “excellence” will be represented quantitatively, 
later on, by the “performance index (PI).” 

Various adjustments of process settings (setpoints, biases, 
manually set, etc.) result in various operating performances 
with the current facilities; optimal adjustment decisions are 
those that yield the most valuable combination of operating 
performance metrics possible. 

Finding and applying these optimal adjustments is some- 
times called operation optimization, process optimization, 
optimizing control, and tuning or fine-tuning the process 
(here tuning the supervisory control, not the regulatory con- 
trol). Optimization — doing as well as possible within the 
scope of decisions and factors of interest — truly makes sense 
only in quantitative terms. The quantitative technologies 
available, each with several alternatives, are first-principle 
models, design of experiments (DOE), sequential empirical 
optimization (SEO), and neural networks (NN). All technol- 
ogies except SEO require initiation by creating a mathemati- 
cal model of the process, which of course will need to be 
updated as the process and business change. 
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This chapter focuses on operation optimization with the 
SEO intelligent analytics applied as macrosupervisory con- 
trol, working above the existing controls, even those with 
various levels of advanced process control. This chapter and 
the references explain why SEO is faster, more comprehen- 
sive, and safer than all other technologies mentioned above. 
In particular, SEO implementations are likely to manage all 
aspects of performance automatically, while the other solu- 
tions are highly complex that they are applied only to a few 
selected performance metrics. 

The area and benefits of production operation optimiza- 
tion (with existing assets) are a vast reservoir of potentially 
enhanced reliability and financial gains at this time. The 
SEO experience indicates that a large fraction of all indus- 
trial production operations underutilize the existing assets in 
terms of the performance metrics selected by the user. The 
benefits achieved with SEO usually range in economic value, 
for high-volume production, from a few $100 K/year per pro- 
cess to a few $1000 K/year per process. 

Operation Optimization (Process Optimization) 

“Regulatory control” is to maintain controlled process vari- 
ables close to their setpoints, where the setpoints are the 
adjusted desired values. Regulatory control is one of the 
functions that can be optimized, with its tools; it contributes 
to better operations, but for the scope of this chapter, it is con- 
sidered different from “operation optimization” (even though 
SEO can be applied to regulatory control as well, especially 
for multiple interacting PID loops). 

The scope of operation optimization (process optimiza- 
tion) is as follows: 

• Given objectives, the current equipment and the 
regulatory control system, and — if relevant — given 
the current values of uncontrolled inputs to the pro- 
cess (imposed conditions such as fuel quality, load/ 
demand, and state of the equipment) 

• Determine the best adjustments of the machine set- 
tings to achieve the most desirable combination of per- 
formance metrics 

In order to perform the above scope of “operation optimi- 
zation” well and practically, when there are more than two 
inputs and more than two outputs, one needs mathematical 
multidimensional modeling and optimization technologies. 

The adjustment of setpoints is a function in supervisory 
control system. If there is already some control logic that 
makes automatic setpoint adjustments depending on condi- 
tions, then the optimization of those adjustments is achieved 
by adding a bias to each control logic setpoint value. 

Note that this scope of operation optimization is different 
but complements the typical engineering improvements by 
changing the equipment, instruments, sensors, and the con- 
trol system. 


SEQUENTIAL EMPIRICAL OPTIMIZATION 

SEO is based on operating data. Therefore, the technology 
will reflect how the process actually operates with all its 
intricacies — as are DOE and NN — however, it is subject to 
limitations in the quality of available data. On the other hand, 
solutions based on first-principle models are subject to limi- 
tations in representing reality sufficiently accurately. 

The principal strength of SEO (as compared to DOE and 
NN) is an immediate use of the information in the sequen- 
tially collected operating data, taking advantage of that 
information as soon as it becomes available to decide how 
next set adjust machine settings. This has some similarities 
to the feedback mechanism in regulatory control, but as we 
shall see below, it is quite different in nature (it is feed-for- 
ward). Because of learning immediately from sequentially 
available information, SEO is also able to safely and gradu- 
ally extrapolate into areas without existing experience. 

Nature of SEO Cycles 

SEO intrinsically solves a multivariable (inputs and outputs) 
problem, and thus all the data needs to be analyzed — and 
displayed — together. Figure 27.1 is a general schematic of the 
SEO cycles, with possible examples of specific variables in a 
variety of industries. The examples at the end of this chapter 
have the variables defined either in an input/output decision 
diagram or in a list of variables. 

SEO operates as follows: starting at the current (base- 
line) process adjustments, it modulates the adjustments (i.e., 
by small changes) and, from the consequent performance 
measurements, learns how to modulate the adjustments next, 
to gradually and sequentially extrapolate — and eventually 
interpolate — and converge to optimal operations. 

In addition, there is an important compounding — 
bootstrapping — benefit to SEO: since it quickly starts con- 
verging to better production runs (as compared to those made 
at the equivalent time with the other methods, some being 
experiments), it also produces useful data closer to the opti- 
mum to learn faster where the optimal combination of adjust- 
ments is. 

It is worth noting that SEO does not require process mod- 
els or historical data to start with. By comparison, all the 
other analytical solutions start by creating process models, 
which involves significant expertise, costs, and time. The 
data necessary for DOE and NN at all likely combination 
of inputs, and equivalent for model validation including for 
first-principle models, are obtained at the price of poor per- 
formance, violating basic requirements, and risks to damage 
the process. To identify the (multiple, as conditions change) 
optimal operations, those models need to be validated where 
the optimal running states are, which is a stringent require- 
ment before doing the optimization analysis — else will need 
to collect more operating data until covering the range of 
optimal adjustments. 
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FIG. 27.1 

SEO cycles diagram. 

Historically, the principles of sequential statistics were 
written by Wald [1], and developed as far as EVOP [2] (evo- 
lutionary operations, a sequential DOE scheme). Separately 
and at about the same time, a heuristic approach (not based 
on analytical model building) was developed called Simplex 
[3], Both the Ultramax SEO described in this chapter and 
the Grabitech’s Simplex are commercial products. Other 
researchers have applied a variety of sequential analysis 
algorithms aimed at quick optimization of various objec- 
tives (e.g.. Ref. [4]) and a list of publications at the end of 
the references. 

The periodic, repetitive steps in an SEO cycle are 

(a) Modulate (multiple) adjustments 

( b) Get (multiple) operating outcome measures 

(c) Enter Run Data into SEO software 

(d) Get from the SEO analysis alerts about the process 
and optimization, and advice for the adjustment mod- 
ulation for the next cycle, together with expected out- 
put values 

The adjustments for the first cycle — the baseline — are the 
previous standard operating procedure process adjustments, 
from which operating improvements are measured. Then 
SEO cycles gradually and safely extrapolate adjustments 
toward the optimum, while all are engaged in regular pro- 
duction, asymptotically approaching optimal performance, 
learning from each new run data how to do better for the 
following cycles. The SEO cycles usually last 30 min to 2 h in 
continuous processes, or the batch time. 


In more detail, what happens in step (d) is as follows: 

• Updates: With the new run data it updates a database 
of historical run data that includes adjusted inputs, 
uncontrolled inputs, and their outcomes or outputs, 
and other variables 

• Learns: With the updated database, it creates/updates 
prediction models for each output as a function of all 
inputs (the models are temporary until the next update) 

• Optimizes: With the updated models, a search engine 
finds the best adjustment modulations for the next 
cycle, the advice 

These three steps also generate alerts about the data, state 
of operations, and state of optimization. These behind-the- 
scenes activities are explained in Section Mathematical 
Structure of SEO. 

A thorough and mature SEO, through automatic data 
exchange data with the process control system, can reli- 
ably run closed-loop — until some anomaly in operations is 
detected, in which closed-loop turns itself off to bring the 
matter to the attention of plant personnel. 

MATHEMATICAL STRUCTURE OF SEO 

A process is conceptually represented by the list of variables 
and their roles (some with lower constraint LC, upper UC, or 
both constraints), one variable being the PI (objective func- 
tion) to be maximized or minimized (after reliability meeting 
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constraints), known calculations, and the sensitivities mini- 
mum important differences (MID, explained below in the 
context of their use). Putting all together results in what is 
called the optimization plan (OP) for operations. 

We will follow the approach of defining, for the sake of 
clarity, what is the meaning of “optimization” with perfect 
knowledge and then addressing a solution to resolve the lack 
of knowledge, which in this Chapter through SEO. 

The following mathematical formulation of opera- 
tion optimization of an I/O process is from the Blue Book, 
Mathematical Formulation of SEO by Moreno [5] (in this 
formulation, lower case bold letters are vectors and upper- 
case bold letters are matrices): 

Ideal Solution with Perfect Knowledge 

Let an operating process behave like 

y = f(x,u) + e (27.1) 

where 

x is the vector of adjusted inputs for each run data set, the 
decisions that affect operational performance 
u is the vector of uncontrolled inputs for each run data 
set that affects operational performance, with values 
determined elsewhere (imposed by people or nature). 
In static conditions, u does not exist 
y is the vector of outputs for each run data set (after any 
transients due to readjustments) that represents the out- 
comes or consequences from the inputs x,u. y 0 is the 
value of the PI or objective function 
f(x,U) is the output mean vector, or the “response.” In con- 
tinuous processes, it is the steady-state after transients 
created by adjustment modulations. The output vector 
f has the indices 0, 1, 2, 3, ... where f 0 (x,u) is the value 
of the PI (objective function) 

f as usual in production processes, is relatively “smooth” 
in the area of interest, f could possibly be changing 
slowly with time (or equivalently, changing because 
of slow changes in unknown uncontrolled inputs not 
included in u). Fast changes are represented by noise 
e is the noise vector produced by a bell-shaped distribution 
with mean zero and covariance matrix 2. The square 
root of the diagonal elements are the standard devia- 
tions, sigma or noise vector n, most likely affected 
by x,u. e is a property of the process and the OP, in 
particular, of the inputs included (e.g., the uncontrolled 
inputs u included) and how y is measured (e.g., aver- 
ages of more raw data may have lower noise, especially 
if u does not change too quickly) 

For this formulation, the concept of noise is smaller than 
that in most quality control analyses because major changes 
in the known uncontrolled inputs u do not contribute to noise, 
while they do in regular methods. 


Note: Fet us understand, contrary to remarks in vogue in the 
last few years in certain fields, that almost no processes pro- 
duce any data following a normal distribution (except aver- 
ages of large amounts of data obtained similarly — due to the 
central limit theorem). There are almost no physical principles 
that indicate that process outputs should be normally distrib- 
uted — an exception being the energy emitted by a black object 
versus the log of the frequency. On the contrary, for instance, a 
process often places limits, resulting in trimming tails, which 
destroys a normal distribution. This has been confirmed by 
experience, where with sufficient data one usually can prove 
statistically that process output distributions are not normal. 

Note that in this model, “errors” are assigned only to the 
outputs. The inputs are presumed to be perfectly correct. 

With this formulation, the definition of optimal opera- 
tions is as follows: 

Maximize the performance index ( the objective function ) — 
while the data satisfies reliably constraint requirements — by 
adjusting the settings x for given values of u. (It could be 
“minimize" if the objective function is a cost or loss function.) 
The constraint requirements are as follows: 

(a) The mean process values seldom violate the 
constraints. 

Note: In this chapter, “mean” is used to express the 
average for the same values of the inputs for one run; 
“average” is used for the average of several runs across 
time (with possible changes in inputs). 

(b) The actual output data for each run is seldom beyond 
a constraint by more than the “minimum important 
difference” (MID, defined by the user). This is the 
mathematical representation of “reliability” meeting 
basic requirements, as the evaluation is made by the 
data, not by the unknown reality. MID can be zero, 
while MID > 0 allows for a commonly understood 
“gray” area given that the “data” are a noisy repre- 
sentation of reality. 

Mathematically, given the value of u, 

Max x/u (f 0 (x,u)} f 0 is the objective function 

s.t. (subject to; i.e., while satisfying these constraint 

requirements) 

The subscript “i” below applies to those variables with 
respective upper or lower constraints 
Requirement (a) 

Xj < UC; (as a decision, it easily obeys constraints) 

X; > LQ 

fi(x,u) < UC; most of the time 
f ; (x,u) > FCj most of the time 
Requirement (b) (for data values) 

y; = f;(x,u) + 8j < UQ + MID; most of the time 
y; = f;(x,u) + 8j > LC; - MID, most of the time 

The “most of the time” requirement is translated below into 
requiring that the condition be satisfied for all but extreme 
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values of le ; l, namely, by le ; l < 3n, (3 -noise-sigmas), called 
the 3-sigma protection (we are not using the criterion of 
probability of violating constraints to avoid having to make 
assumptions about the distribution of the data). So, 
y i =f i (x,u) + 3n, < UC ; + MID; 
y i =f i (x,u) - 3nj > LQ - MID; 

So, bringing both constraint requirements (a) and (b) into 
one composite set of equations, the optimum is defined as 

Max x /u {f 0 (x,u)i 
S.t. X; < UC; 

X; > LCj 

fj(x,u) < UC; - max{3n; - MIDj, 0} 
fj(x,u) > LQ + max(3n ; - MID ; , 0) 

which defines the desired optimum adjustments x* and the 
mean optimum outputs y* = f(x*,u); or [x*,y*] for each value 
of u. 

The last two inequalities define the effective operat- 
ing area for the mean outputs fi(x,u), where the safety 
buffer = max{3n i -MID i , 0} is illustrated in Figure 27.2. If 
the safety buffer is larger than zero for active constraints, 
this constitutes a loss in operating area that results in a loss 
of opportunities for higher f 0 . 

The practical objective is not just to find the optimal x 
but also to define the set of xs that satisfy all constraints and 
produces mean outputs no further away from the optimal 
f 0 (x*,u) than MID 0 . The set of such xs defines the window of 
operations or the “sweet operating area.” 

If there are outputs with upper and lower constraints, a 
sufficiently large safety buffer could result in, there being no 
operating range, and then the process is totally inefficient. 



FIG. 27.2 

The “effective operating area” to satisfy (b) requirement; the area 
of acceptable output means. 


In this case, optimization would just reduce constraint viola- 
tions, and the problem needs to be fixed, probably with capi- 
tal expenses involved, such as the first point below. Following 
are other procedures to reduce losses in operating area due 
to high noise: 

• Find the cause of noise and reduce or eliminate its 
effects (reduce noise), such as through better regula- 
tory control and/or more accurate sensors 

• Find the cause of the noise and include it as a mea- 
sured uncontrolled input (recognize the source of its 
effects in the analysis, eliminating it from the noise, 
and enable SEO to compensate for its known value) 

• Relax requirements: relax the constraints and/or 
enlarge the MIDs 

Sequential Empirical Optimization 

This section deals with how to approach the ideal solution 
described above when there is a lack of knowledge about 
both the form and coefficients of f(x,u) and of E — except for 
knowing that since it is a production process, the response is 
reasonably “smooth.” 

As explained above, the SEO cycle basically 

• Adjusts the production process with x and gets one run 
of operating data {x,u,y} t 

• Enters the data set into the software and gets in return 
a safe x for the next run 

The sequential analysis above emulates the process of a mind 
gaining experience through repeated action, and using the 
remembered actions, conditions, and outcomes to make bet- 
ter action decisions in the future — an aspect of learning and 
skill development. 

Sequential analysis has an advantage [1] that it uses the 
information in newly collected run data immediately in order 
to refine knowledge and to decide where to run next, thus 
increasing value generated right away. By comparison, in tra- 
ditional more parallel empirical studies (e.g., DOE, NNs), the 
information is not obtained until the data are analyzed at the 
end and turned into increased value even later. The simple 
advantage of early use of information is the most important 
reason why a properly implemented SEO technology is the 
fastest empirical optimization available today. In addition, 
SEO is faster due to bootstrapping, that is, getting data close 
to optimal sooner helps create better prediction models as to 
where the optimum is. 

The gradual SEO approach used during regular produc- 
tion could converge to a local optimum, as it would be unde- 
sirable to go searching down a “valley” of lower operating 
performance hoping to find a larger “mountain” on the other 
side. Finding such other separate mountains to climb — if 
they exist at all — belongs to engineering, process develop- 
ment, or R&D; it should not be attempted as part of daily 
production. 
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The performance of a SEO solution is evaluated by 

• How quickly the series {x,u,yj t converges to the opti- 
mum {x*,u,y*} 

• How closely the series {x,u,y} t converges the optimum 
{x*,u,y*} 

• In production, which is a case of continual value gen- 
eration, the evaluation is the cumulative value gener- 
ated since the starting of optimization project 

For instance, the cumulative value generated with SEO will 
be compared with the poor performance (and process abuse) 
while collecting experimental data for technologies such 
as DOE and NNs and the ongoing old performance of runs 
with baseline adjustments when the first principle models 
are designed and implemented. In addition, for all of these 
approaches, there will be runs with widely distributed adjust- 
ments to validate the models created. 

Depending on the SEO technology used (or on how it is 
adjusted), achievements in terms of the objective function f 0 may 
frequently be increased at the sacrifice of increasing the inci- 
dence of violating some constraints, and thus the SEO needs to 
proceed carefully. The sequence of x t needs to be “intelligent”: 
low risk, effective, and responsive to changes in conditions. 

SEO Analytics 

A SEO solution is based on creating prediction models, as 
described in Section Nature of SEO Cycles: every time a new 
adjustment advice is required: 

• Create/update prediction model(s) m(x,u), which is a 
suitable approximation to f(x,u), based on the avail- 
able data in the database 

• Generating adjustment decisions (advice) by optimiz- 
ing the prediction models 

Note: a technique that is sequential but is not based on predic- 
tion models is Simplex [3], In this chapter, we continue with 
creating model-based predictions that are necessary to con- 
trol the reliability/risk of future data violating constraints and 
necessary to enable dynamic optimization, where adjustments 
compensate for changing values of the uncontrolled inputs. 

A method to create m(x,u)sf(x,u) is least-sum-of- 
squared-errors regression analysis based on the set of accu- 
mulated historical run data {x,u,yj t . The m is either linear or 
quadratic, depending on the amount of data and the quality 
of the fit. At the same time as the models m are fitted to the 
data, the SEO technology also creates 

• An area of confidence (AOC), which is the region in 
x,u where m(x,u) is most accurate — within the region 
covered by the sets {x,u} t . The AOC can be concave 
or even composed of disjointed areas, which can be 
resolved by the proper use of the Mahalanobis dis- 
tance. Optimum estimates and adjustment advices 
are given only within the AOC — where an optimum 


estimate at the edge of the AOC is an indicator that 
optimal operating conditions are seen outside the 
AOC; thus, optimization involves extrapolation. 

• An estimate of the covariance matrix X of e, and 
therefore of the noise vector n. 

Let us recognize two properties about the required model m 
for SEO: 

1 . It is not required to understand the effects of each input 
separately; it is just necessary to be able to predict 
results in the AOC region. Thus, with the protection 
afforded by the AOC, confounding of input effects is 
of little consequence. 

2. Trying to make m accurate away from the optimum 
(e.g., by making its generic form too involved or by 
fitting data away from the optimum) is not only unnec- 
essary but also distorts the fitted models resulting in 
lesser accuracy around the optimum. 

These two properties present less demands on m(x,u) than 
the usual requirements for DOE (especially orthogonality), 
NNs, and first principle models (each for different reasons). 
Thus, when not using user-supplied calculations, the generic 
mathematical models of m are at most a quadratic approxi- 
mation (second-order Taylor) of the response surface f(x,u), 
the “quadratic” being necessary to represent interactions 
(nonlinearities) between inputs. 

If using regular regression to create the prediction models 
m, this leaves the creation of models unresolved in the first 
N + lew cycles, where N is the number of inputs, in which 
case, it will create a singular matrix (X'X) impossible to 
invert. Thus, one solution is to follow the general approach of 
DOE and NN and start with a set of experimental runs, and 
then follow through with cycles; and indeed this has been 
done to refine the results of studies already started with DOE 
or NN. See below for the Bayesian regression alternative. 

The method to create the advice x for the next run based 
on m, n, and AOC is the solution to 

Max x/u {m 0 (x,u)} 0 for objective function 
s.t. Xj < UQ i = inputs with upper constraint UC 

X; > LC; i = inputs with lower constraint LC 

m i (x,u)<UC i -max{3*n i -MID l , 0} i = outputs with 

upper constraint UC 
m i (x,u)>LC i +max{3*n i -MID i , 0} i = outputs with 

lower constraint LC 

{x,u} e AOC 

which is done with nonlinear programming, with this 
prioritization: 

First: inputs: obey their constraints 

Second: when giving advice only (not for the optimum esti- 
mate), give advice where each inputs adjustment modulation 
(change from last value) is not larger than its MID 
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Third: {X,U} e AOC, to assure relative certainty of estimates 

Fourth: outputs obey their constraints (note, the PI is ignored) 

Fifth: optimize the PI (objective function) and sometimes 
add a perturbation to generate data with more information 
for future model building 

Some specific technical features of SEO software (e.g., 
in Ultramax®): 

• Dynamic optimization (as indicated by the presence of 
u in the math models above): Adjustments can depend 
on the current values of important uncontrolled inputs, 
which are the conditions imposed on production. Also, 
it responds quickly to changes in economic factors 
through the recalculation of the PI with new factors. 

• The method includes a procedure to evaluate and 
determine the balance and trade-offs of several per- 
formance metrics (to create the PI) that effectively 
reflects the culture and values of plant personnel. 
When applicable (technically and motivationally), the 
balanced evaluation is determined by the impact of 
the P&L statement by the scope of this optimization, 
that is, letting profit maximization drive directly the 
optimization. 

• To kick-start getting improvements with little accumu- 
lated data generates useful prediction models by using 
Bayesian regression (while other empirical technolo- 
gies would be at the stage of running experimental 
runs). For the mathematical solution, see Colosimo 
and del Castillo [6], 

• The users are not compelled to follow the advice for 
new adjustments, they can modify it, thus keeping the 
user in full control if they choose to, another aspect 
of safety especially in the early stages with little data 
when SEO did not accumulate sufficient amount of 
knowledge yet. 

• The software accepts outputs that are calculated (as a 
function of other variables). Along the same lines, it 
accepts prior process models, even though incomplete 
or inaccurate, useful for faster scaling up. 

• To create more accurate models when there is a lot 
of data, it creates locally accurate models through 
weighted regression — thus avoiding the distortions 
of fitting data away from the region of adjustments of 
interest. 

• This SEO software can be used both with manual data 
entry and with automated data exchange with the pro- 
cess control system, such as a distributed control system 
(DCS) or historian, for up to autonomous closed-loop 
dynamic optimization. This automatic handling of all 
inputs and outputs is the reason why each application 
can be so easily as comprehensive as the user wishes. 

• Since the SEO macro-supervisory-control installation 
does not require a change to the process or controls, it 
can always be “unplugged”; the user is not committed 
to its continued use. 


With this version of SEO, only about 1% of the applications did 
not deliver improvements because, to start with, baseline opera- 
tions were already close to optimal. The rest got improvements 
of the PI, reduction of risks violating constraints, or reducing 
the level of constraint violations, (see examples below.) 

APPLICATION OF SEO 
Procedure 

• A team (or an individual) formulates the OP: 

• Defines objectives for improvements together with 
their performance metrics 

• Defines minimum requirements (metrics and their 
constraints) — such as safety, quality specifications, 
regulations, commitments such as budget, and so 
on, which are the basic and highest priorities to 
be met. 

• Define the process inputs likely to affect the per- 
formance metrics: the adjustable to adjust (manu- 
ally or at the control panel), the adjustable to keep 
constant, and the uncontrolled (imposed conditions 
that one has to work with). Indicate constraints rep- 
resenting acceptable values, as applicable. 

• For when the constraints are satisfied reliably, define 
a PI to be maximized or minimized, which could 
be a single output or a calculated overall operating 
performance by the balance of several performance 
metrics. An ideal PI is, when possible to calculate, 
the total impact economic on the P&L statement. 

• Formulating the OP for production process, docu- 
menting it, and getting management approval takes 
about 2 days; for automatic data exchange, add 
about a week for installing and testing an interface 
of SEO with the control system. 

• Engage in production operations managed with SEO 
cycles as described previously. Follow advice and 
respond to alerts (such as multivariable statistical pro- 
cess control). 

• Periodically document progress. 

• When process or objectives change, update the OP. 

• Do sensitivity analysis of constraints to elicit ideas for 
equipment improvements. 

A detailed description of procedures and specifications are 
explained by Moreno [5]. 

Example Applications 

The six examples below include a description of the process 
and the main results obtained. 

Recall that SEO is a multivariable task, one needs to be 
aware of all the variables involved — easy for the computer, 
more difficult for us — and of course, about the changes in 
all the variables (as improvement in one variable would 
be meaningless if other variables diminish performance 
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significantly). Some ways to convey the multidimensional 
reality are as follows, illustrated in various examples below: 

• See all the variables in a decision I/O diagram. 

• In the plots, the color of each dot represents the worst 
constraint status among all the variables. A red dot 
means that the variable with the worst constraint sta- 
tus had an important constraint violation, orange that 
it had an unimportant violation (<MID), yellow it sat- 
isfies constraints but is close to it, green that all vari- 
ables are safely away from the constraints, and black 
that it is in the “sweet area of operations.” 

• There is a comparison report that looks only at the data 
in one set of runs versus another. The last example in 
Figure 27.10 shows the comparison of last 20 runs under 
the “regime” of SEO with the first run at baseline adjust- 
ments. The right column has a statistical analysis indicat- 
ing the certainty of the changes in average performance 
as supported by the data. Values changes above 3 -sigma 
are considered “definite change” (with very high statisti- 
cal confidence), and the confidence is gradually reduced 
to “undetermined” for values below 1-sigma. 

In the plots, each dot represents the results from one set 
of adjustment modulation. Often the value plotted is the 
10-moving-average, to see trends more clearly. 

Example 27.1: Food Processing — 
Continuous Cooking Ovens 

Description: The process uses a series of steam and gas 
fired stages in a continuous oven to cook unbreaded pat- 
ties. Critical control points and quality are measured by 
internal cooked temperature and weight, as well as color 
of the patty (against a scale represented by pictures). 


Situation: Capacity constrained production — can sell 
more if it is produced. USDA has regulatory specifica- 
tions. The scope of operation optimization is illustrated 
in Figure 27.3. 

Objectives: Maximize profitability impact rate (within 
the scope of operation optimization), with reliability ful- 
filling all requirements. Also, instill in plant personnel a 
quantitative awareness of production. 

Results summary: Economic gains (Tot Sav) were about 
SlOlO/h of production (with conclusive statistical confi- 
dence at 3.5 sigmas), mostly through production rate (and 
Sales) increasing 40% (while satisfying all requirement 
constraints such as critical control points and product 
quality). Yield increased by 0.4% but with no statistical 
confidence yet. The rate of improvements with SEO is 
displayed in the plots of Figure 27.4. 

Except for the first plot, all others are 10-moving- 
average, to perceive changes more readily. 

Example 27.2: Paper Industry — 
Energy Saving 

Process: Tissue paper mill, with through airflow and 
Yankee drying. The scope of operation optimization is 
illustrated in Figure 27.5. 

Objectives: Reduce energy consumption costs while main- 
taining quality constraints, for constant running speed. 

Results summary: First stage: 7.1% reduction in energy 
cost. The optimum was at wide-open recycled airflow. 
Software sensitivity analysis, confirmed by engineering, 
resulted in the desirability to increase recycle fan and duct 
sizes ($50 K). The engineering change was implemented. 
The next optimization, not at the max recycle capacity, 
saved another 8.5% energy cost. Total savings of 15.6% 
were obtained. The findings were applied to 11 other 
mills (total gains were a few million dollars per year). 


Decision input/output diagram 


Adjusted inputs 


Outputs. Outcomes 

5 Steam rates (lb/h) 

Process: 

Patty temp (avg and var) ^ 

> 

> 

2 Oven temps (°F) . 


Patty color (scale) 

> 



Formax speed (spm) 


Patty weight (lb/8 pat) ^ 

— ^ 


Yield (%wt) v 

Uncontrolled inputs 


>> 



Production (klb/h) ^ 

Heat recv flow (ldb/h) ^ 


>> 

> 


Economic variances (K$/h) 

Raw Patty temp. (°F) ^ 

Run time: 60 min 

>> 


Analysis: 15 min 

Perf. Index: Tot savings (K$/h) 

Raw Patty Wt. (lb/patty) ^ 


> 


FIG. 27.3 

Decision input/output diagram, part of the OP. 
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Decision input/output diagram 


Adjusted (control) inputs 


Back press, damper 

— > 

Make-up air damper 


Ratio of two biases 


Exhaust damper 

r 

Uncontrolled inputs 


Furnish characteristics 

— > 


Process: 


Tissue paper 
mill 


Run time: 50 min 
Analysis: 25 min 


Outputs. Outcomes 

Moisture 

Temperatures 

Gas consumption 
Electric power consumption 

Recycled air flow 

Tensile strength 

Caliper 

Score: Energy cost, min 


> 

> 

> 

> 

> 

> 

> 

> 


FIG. 27.5 

Paper energy savings, decision input/output diagram. 

Example 27.3: Plastic Forming — 
Injection Molding/Large Parts 

Process: The operation, set up in the machine maker's 
test facility, used a very large twin-screw injection mold- 
ing machine. It was equipped with a flash gate mold for 
making autofenders. Quality was defined by a combina- 
tion of surface and dimensional requirements. The scope 
of operation optimization is illustrated in Figure 27.6. 

Objectives: Minimize deviation from dimensional targets 
and improve surface quality index. 

Results summary: Superior quality parts were produced 
after only a few days of testing enabling the start of 


production runs based on the developed procedure. Also 
further analysis of the data for optimal adjustments of the 
models for the gap test points disclosed that the top cor- 
ner and pillar required three times the “hold time” of the 
optimal overall part. This prompted the mold makers to 
open the flash gate at these points to cause greater plastic 
flow to these areas of the mold. 

Example 27.4: Power Generation, Lignite 

Situation: Litvinov Boiler K19 is one of four new boilers, 
500klb/h, using Lignite, ABB low-NO, corner burners, 
with Honeywell’s TotalPlant® Solutions System plus the 


Adjusted (control) inputs 
Injection speed (s) 

Injection pressure (PSI) 
Hold pressure (PSI) 

Shot size (in.) 

Melt temp. (°F) 

Mold temp. (°F) 

Hold time (s) 

Cure time (s) 


Decision input/output diagram 


Process: 


Injection 
molding/ 
large parts 


Run time: 24 h 


Outputs. Outcomes 
Part weight (kg) 


Shrinkage (%) 


Surface quality knit 


Surface quality lifter 


Dimension 1 (gap) 


Dim 5,7,9, flipr, bmpr 


Perf. index: Overall reduction of 


variations from targets 


> 

> 

> 

> 

> 

> 

> 


FIG. 27.6 

Paper energy savings, decision input/output diagram, part of the OP. 
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Control inputs 
0 2 flue gas 
Tert ab split 
Tert lowR 
Tert total R 
Sec AB split 
Mill Bias 1 
Mill Bias 2 
T Mill out 
Uncontrolled inputs 
Steam prod 
PAFR 


Chemopetrol 

Litvinov 

Boilers 17-20 


Boiler efficiency 
improvement 


Outputs 

Efficiency 


T flue gas 


CO 


NO v 


T Mill 1 


T Mill 2 


T Mill 3 


F spray 


Feeder 1 


Feeder 2 


Feeder 3 


FIG. 27.7 

Chemopetrol plant, decision input/output diagram, part of the OP. 


“Individual Boiler Optimization®” (IBO), which includes 
this SEO. IBO provides a real-time shell providing auto- 
matic operations management either open (Advisory) or 
closed-loop. For more details, see www.ultramax.com/ 
WhitePapers/Closed-loop_Optimization.pdf. The scope 
of operation optimization is illustrated in Figure 27.7. 

Objectives: Maximize efficiency while control NO* 
within regulations and satisfy all safety constraints. 

Results: The table and plots show the first 160 cycles that 
took about 120h of operations. Efficiency increased by 
4.3%; NO r and CO increased but were comfortably below 
the maximum value allowed (CO<250mg/Nm 3 and 
N0,<500mg/Nm 3 ). The rate of improvements with SEO 
is displayed in the plots of Figure 27.8. 

At the beginning (in the plots, not the first run at 
baseline adjustments in the table), there were violations 
of input constraints generated by plant personnel overrid- 
ing SEO's advice (which never violates input constraints), 
and after run about run #100, the SEO advice was fol- 
lowed well, running closed-loop. This is illustrated in 
Mill Bias 1 plot (note that it is has the individual cycle 
(run) values rather than the moving average). The same 
happened with Mill Bias 2. In any case, these adjusted 
inputs did not change much from the baseline (the first, 
reference) run. To be consistent, 0 2 flue gas is also plotted 
with individual cycle values. 


Results: All constraints were satisfied. The table and 
plots (10-moving-average) show the first 665 cycles that 
took about 400 h of operations. Efficiency increased 
by 1.5% (as uncontrolled inputs coal heat increased 
by 11% and Coal H 2 0 increased by 8.1%); LOI was 
reduced by 49%. To see the improvements in efficiency 
above the increase in coal heat, in the second plot in 
Figure 27.9 “efficiency versus coal heat,” the efficiency 
behaved in a particular manner in the first -180 runs, 
and then operations gradually shifted to a higher effi- 
ciency pattern. 

Example 27.6: Public Services 
Electric and Gas — Selective 
Noncatalytic Reduction 

The project included combustion optimization of boiler 
Mercer Unit 2 and performance optimization of the selec- 
tive noncatalytic reduction (SNCR) unit that treats Unit 
2’s exhaust gases running at full load. Phase 1 optimiza- 
tion of Unit 2 resulted in reductions of heat rate by 1.2% 
and NO, by 10%. The following is the optimization of the 
SNCR in phase 2. 

The SNCR is a NALCO unit. Both the boiler and the 
SNCR units have their own DCS. The SNCR unit focused 
on zones 5, 6, and 8. The primary adjustments were made 
to air and water pressure and urea volume. SNCR com- 
parison table of last 20 runs versus first run at baseline is 
given in Figure 27.10. 

Objectives: Reduce ammonia slip (economic con- 
siderations) while satisfying safety and performance 
requirements. 

Results: All constraints were satisfied. Ammonia was 
reduced 80%, from 7.1 to 1.4ppm. NCT, and opacity were 
reduced by statistically inconclusive amounts. CO was 
reduced to an almost certain amount of 60%. The rate 
of improvements with SEO is displayed in the plots of 
Figure 27.11. 

The large modulation (readjustments) from the first to 
the second cycle did override SEO’s advice, which would 
not have moved so fast. This reflects awareness on the 
part of plant personnel that could not be implemented 
before, and is an example of the feature in SEO that lets 
the users make the readjustments they wish and SEO still 
learn from the consequent results. It is also an indication 
that SEO is a toolset to expand the capabilities of the user, 
not to substitute the users. 


Example 27.5: Power Generation, Coal 

Situation: Zou-Xian Unit #6 is a 550 MW Dong Fang, 
tangential, coal fired boiler. It has five mills, two sets of 
three levels of burners, and two sets of four levels of sec- 
ondary air feeders, controlled by a Westinghouse WDPF 
DCS. The quality of coal varied considerably. 

Objectives: Given safety/reliability considerations, the 
main performance objective was to maximize efficiency 
and reduce LOI (for quality of ash). 


CONCLUSIONS 

SEO is a most effective (fast, comprehensive, and safe) 
toolset to achieve improvements and optimization in opera- 
tions with the existing processes, as evaluated according to 
the user’s metrics, with data provided by the user. There is 
an established track record in several industries including 
chemical (continuous and batch), power generation, phar- 
maceuticals, paper and nonwovens, food processing, and 
automotive. 
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Run number— Efficiency run data Run number— CO run data 



0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 

Run number Run number 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — u.uz, 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 


Run number Run number 


These two plots are individual runs, rather than moving average 
Run number— 02_Flue„gas run data Run number— Mill_bias_l run data 



0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150 160 


Run number Run number 

FIG. 27.8 

Chemopetrol plant, rate ofSEO improvements. 
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COMPARISON Report/CHANGE Analysis 

APPLICATION: C : \UMC\ ... Mersncr3 . Umax | | Mercer SNCR Test with 5, 6 & 8 
ULTRAMAX® 7.1.h ©1982-2010, Ultramax Corporation. All rights reserved. 


Filters for FIRST (REFERENCE) SET : Run range: 1 1 

Filters for SECOND SET : Run range: 53 72 

First Set Second Set DIFFERENCE 
| AVG n | AVG n | Delta % Stat . Signif ic 
III (# Sigmas) 

CONSTRAINT VIOLATIONS: 


% runs w/some INPUT viol ! 
% runs w/some OUTPUT viol ! 


All input constraints satisfied. 
All output constraints satisfied. 


VARIABLE VALUES (averages) : 


# NAME UNITS ROLE 

Adjusted (decision) inputs . 

First 

1 

Second 

1 

Delta 

% 

Sigmas 

1 

Zone 6 wat 

psi 


1 

60.0 

1| 

74.9 

20 | 

14.9 

25.% 

3.0 

2 

Zone 6 ure 

Urea 

gpm 

1 

99.00 

1| 

95.80 

20 | 

-3.20 

-3.2% 

-1.4 

3 

Zone 5 wat 

psi 


1 

80.00 

1| 

84.15 

20 | 

4.15 

5.2% 

1.8 

4 

Zone 5 ure 

Urea 

gpm 

1 

99.00 

1| 

100.40 

20 | 

1.40 

1.4% 

0.8 

5 

Zone 8 wat 

psi 


1 

65.00 

1| 

73.60 

20 | 

8.60 

13.% 

2 . 9 

6 

Zone 8 ure 

Urea 

gpm 

1 

113.0 

1| 

92.2 

20 | 

-20.8 

-18.% 

-3.0 

7 

Zone 6 air 

psi 


1 

40.00 

1| 

42.15 

20 | 

2.15 

5.4% 

1.8 

8 Zone 5 air psi 1 

Performance Index (Minimize) 

40.00 

1| 

38.25 

20 | 

-1.75 

-4.4% 

-0.9 

10 Ammonia ppm 6 

Other consequences , outputs 

7.10 

1| 

2.04 

20 | 

-5.06 

-71.% 

-2.7 

9 

NOx 

ppm 


5 

710.0 

1| 

697.0 

20 | 

-13.0 

-1.8% 

-0.3 

11 

CO 

ppm 


5 

8.30 

1| 

3.43 

20 | 

-4.87 

-59.% 

-3.0 

12 

NOx lbs 

lbs/MBtu 

5 

1.2800 

1| 

1.2740 

20 | 

-.0060 

-0.47% 

-0.1 

13 

02_Stack 

% 


5 

7.0600 

1| 

7.0150 

20 | 

-.0450 

-0.64% 

-0.1 

14 

Opacity 

% 


5 

6.160 

1| 

5.560 

20 | 

-0.600 

-9.7% 

-0.6 

15 

Chem flow 21 

GPH 


5 

314.0 

1| 

286.5 

20 | 

-27.5 

-8.8% 

-2 . 6 

16 

Chem Flow 22 

GPH 


5 

314.0 

1| 

286.5 

20 | 

-27.5 

-8.8% 

-2 . 6 

17 

02 Econ SH 

% 


5 

3.900 

1| 

3.717 

20 | 

-0.183 

-4.7% 

-0.3 

18 

02 Econ RH 

% 


5 

3.800 

1| 

3.908 

20 | 

0.108 

2.8% 

0.3 

19 

Water_SH_5 

gp m 


5 

12.000 

1| 

12.825 

20 | 

0.825 

6.9% 

2.0 

20 

Water_SH_6 

gpm 


5 

9.00 

1| 

11.40 

20 | 

2.40 

27.% 

2.4 

21 

Water_SH_8 

gpm 


5 

8.00 

1| 

10.05 

20 | 

2.05 

26.% 

2.1 

22 

Water_RH_5 

gpm 


5 

13.000 

1| 

13.150 

20 | 

0.150 

1.2% 

0.4 

23 

Water_RH_6 

gpm 


5 

9.00 

1| 

10.25 

20 | 

1.25 

14 .% 

1.0 

24 

Water_RH_8 

gpm 


5 

8.000 

1 1 

8.850 

20| 

0.850 

11.% 

1.0 


FIG. 27.10 

SNCR comparison table of last 20 runs versus first run at baseline. 
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